Activation barrier scaling and crossover for noise-induced switching in a 
micromechanical parametric oscillator 
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We explore fluctuation-induced switching in a parametrically-driven micromechanical torsional 
oscillator. The oscillator possesses one, two or three stable attractors depending on the modulation 
frequency. Noise induces transitions between the coexisting attractors. Near the bifurcation points, 
the activation barriers are found to have a power law dependence on frequency detuning with critical 
exponents that are in agreement with predicted universal scaling relationships. At large detuning, 
we observe a crossover to a difi'erent power law dependence with an exponent that is device specific. 
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Fluctuation-induced escape from a metastable state 
plays an important role in many physical and biologi- 
cal phenomena. As shown in Kramers' early work for 
systems in thermal equilibrium the rate of escape 
depends exponentially on the ratio of an activation en- 
ergy to the temperature. Activated transitions have been 
studied in depth for systems in thermal equilibrium. Re- 
cently, there is much interest for activated transitions in 
systems that are far from equilibrium. A number of sys- 
tems, such as electrons in Penning traps |3, Josephson 
junctions 0], micro and nanomechanical oscillators 0,0] 
and atoms in magneto-optical traps 0, 0| , develop mul- 
tistability under sufficiently strong periodic driving, but 
are monostable otherwise. These systems are not char- 
acterized by free energy and the transition rate must be 
calculated from system dynamics 121; 1^ ^| ■ 

In both equilibrium and non-equilibrium systems, as a 
system parameter /i approaches a bifurcation value /i;,, 
the activation barrier decreases to zero and the number 
of stable states of the system changes. In general, the 
activation barrier is determined by the device parame- 
ters and depends on the specifics of the system under 
study. However, at parameter values close to the bifur- 
cation point, the activation barrier is expected to ex- 
hibit system-independent scaling. The activation barrier 
is proportional to |/i — with a critical exponent ^ that 
is system independent and depends only on the type of 
bifurcation For instance, in a Duffing oscillator reso- 
nantly driven into bistability, spinodal bifurcations occur 
at the boundaries of the bistable region. One stable state 
merges with the unstable state while the other stable 
state remains far away in phase space. Recent exper- 
iments in micromechanical oscillators and rf-driven 
Josephson junctions ^3 have confirmed the theoreti- 
cal prediction j^, ^3 ^^^^ activation barrier scales 
with critical exponent 3/2 near spinodal bifurcations in 
driven systems. On the other hand, a different critical 
exponent of 2 is expected at a pitchfork bifurcation in 
systems where all three states merge Such bifur- 

cation commonly takes place in parametrically driven 
systems where period doubling occurs. For instance. 



fluctuation-induced phase slips were observed in para- 
metrically driven electrons in Penning traps between 
two coexisting attractors and transitions between three 
attractors were studied in modulated magneto-optical 
traps 1^. To our knowledge, the activation barriers have 
not been measured over a wide enough parameter range 
in these parametrically driven systems to demonstrate 
the universal scaling at driving frequencies near the two 
critical points and the crossover to system-specific depen- 
dence at large frequency detuning. 

In this Letter, we report measurements of the ac- 
tivation barrier for fluctuation-induced switching in a 
parametrically-driven micromechanical torsional oscilla- 
tor, a system that is far from thermal equilibrium. The 
spring constant of our device is modulated electrostat- 
ically near twice the natural frequency. Under suffi- 
ciently strong parametric modulation, two pitchfork bi- 
furcation points exist. At the supercritical bifurcation, 
there emerge two stable oscillation states that differ in 
phase by tt. At the subcritical bifurcation, an addi- 
tional, stable state with zero oscillation amplitude ap- 
pears. Noise induces transitions between the coexisting 
attractors. By measuring the rate of random transitions 
as a function of noise intensity, we deduce the activation 
barrier for switching out of each attractor as a function of 
frequency detuning. Near both bifurcation points, the ac- 
tivation barriers are found to depend on frequency detun- 
ing with critical exponent of 2, consistent with the pre- 
dicted universal scaling in parametrically driven systems 
jl4l | . Away from the immediate vicinity of the bifurcation 
point, universal scaling relationships for the activation 
barrier no longer hold. We find that in our parametric 
oscillator, the dependence of the activation barrier on fre- 
quency detuning changes from quadratic to 3/2**^ power. 
As the driving frequency increases further, the activation 
barrier becomes independent of the driving frequency for 
switching between the two oscillation states. In contrast, 
the barrier for switching out of the zero-amplitude state 
increases monotonically with frequency detuning. As a 
result, the occupation of the zero-amplitude state even- 
tually exceeds the two oscillation states. 
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FIG. 1: Response of the oscillator at a;/2 versus the frequency 
of parametric driving u. The solid and dashed lines represent 
the stable attractors and the unstable oscillation states re- 
spectively. Inset: cross-sectional schematic of the device (not 
to scale). 

The micromechanical torsional oscillator in our exper- 
iment consists of a movable polycrystalline silicon plate 
(500/im X 500/im) suspended by two torsional springs, as 
shown in the schematic in the inset of Fig. ^ Two elec- 
trodes are located underneath the top plate. A periodic 
driving voltage Vd = Vdc + Vac cos ut is applied to one of 
the electrodes, where the driving frequency uj = 2uJo + t is 
close to twice the natural frequency Wo(20581.7 rad s~^). 
The top plate is therefore subjected to a periodic elec- 
trostatic torque, the angular gradient of which modulates 
the spring constant. The equation of motion is given by: 

e + 27^ + ^ cosujt^ + a9^ + 139^ = (1) 

where = C" {9o)VdcVac is the amplitude of modulation 
of the spring constant, C"{9o) is the second derivative 
of the capacitance between the driving electrode and the 
top plate with respect to 9 evaluated at the equilibrium 
angular position 7 is the damping constant, / is the 
moment of inertia, a and [3 are the nonlinear coefficients. 
Torsional oscillations of the top plate are detected capac- 
itively by the other electrode. All measurements were 
performed at 77 K and < 10^^ torr. The Q of the os- 
cillator is sufficiently high 7500) so that the response 
of the oscillator at ^ 2wo is negligible. Details of the 
oscillator can be found in Ref. |l5| . 

When the amplitude of the spring modulation exceeds 
a threshold value fc^ = 4^07/, period doubling occurs 
[1^ . Oscillations are induced at half the modulation fre- 
quency in a range close to ojo- As shown in Fig. Q there 
are three ranges of frequencies with different number of 
stable attractors, separated by a supercritical bifurcation 
point Wfci — 2ujo + and a subcritical bifurcation point 
= — (^p, where Wp = \Jk'^ — k'^/2Iujo. In the first 



region (w > wji ~ 41174 rads"""^), no oscillations take 
place, as the only stable attractor is a zero-amplitude 
state. At cjfti, there emerge two stable states of oscil- 
lations at frequency cj/2 that differ in phase by tt but 
are otherwise identical, because of the symmetry with 
respect to a translation in time by 2ti/lo. These two sta- 
ble states are separated in phase space by an unstable 
state with zero oscillation amplitude (dashed line in Fig. 

At frequencies below LLib2{^ 41150 rad s~^), the zero- 
amplitude state becomes stable, resulting in the coexis- 
tence of three stable attractors. These stable states are 
separated in phase space by two unstable states indicated 
by the dashed line in Fig. ^ 

The presence of noise allows the oscillator to occasion- 
ally overcome the activation barrier and switch between 
the different attractors. Since the parametrically driven 
oscillator is far from equilibrium and is not characterized 
by free energy, the transition rate cannot be determined 
from the height of the free energy barrier. Theoretical 
analysis indicates that the transitions remain activated 
in nature 0: 

r cx e^/-f" (2) 

where V is the transition rate and Im is the noise inten- 
sity. In general, the activation barrier R is determined 
by the device parameters such as the damping constant, 
nonlinearity coefficients and the driving frequency. Near 
the bifurcation points, the system dynamics is charac- 
terized by an overdamped soft mode and R decreases to 
zero according to |a; — Wbj^, where the critical exponent ^ 
is universal and depends only on the type of bifurcation. 
In a parametric oscillator, the supercritical and subcrit- 
ical bifurcations involve the merging of two stable oscil- 
lation states and an unstable zero-amplitude state (at 
LO}ji) and the merging of two unstable states and a zero- 
amplitude stable state (at iOb2) respectively. When three 
states merge together in such pitchfork bifurcations, the 
critical exponent is predicted to be 2 Away from 

the bifurcation points, the scaling relationship no longer 
holds and different exponents were obtained depending 
on the nonlinearity and damping of the system. 

In order to investigate the transitions between stable 
states in our parametric oscillator, we inject noise with a 
bandwidth of 600 rad centered at ujo- Figures 2a and 
2c show respectively the oscillation amplitude and phase 
at a driving frequency in the range of two coexisting at- 
tractors. Transitions can be identified when the phase 
slips by TT. The two oscillation states have the same am- 
plitude. These two attractors can also be clearly iden- 
tified in the occupation histogram in Figs. and ^p. 
Figuresl^D and|2Ji show switching events at a driving fre- 
quency with three attractors, where the zero-amplitude 
state has also become stable. In contrast to Fig. Et, 
the oscillator switches between two distinct amplitudes. 
At high amplitude, the phase takes on either one of two 
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FIG. 2: Oscillation amplitude (a) and phase (c) for a driving 
frequency of 41163.0 rad s~^. For driving frequencies between 
ujbi and UL!b2, transitions occur when the phase slips by tt. (b) 
When the driving frequency (41124.8 rads~^) is lower than 
uJb2, transitions involve jumps in the amplitude, (d) In the 
high amplitude state, the oscillation phase takes on either 
one of two values that differ by tt. The phase fluctuates when 
the oscillator is in the zero-amplitude state. 
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FIG. 3: The occupation in phase space at four different ui^s. 
X and Y denote the two quadratures of oscillation. The grey 
scale represents the number of times that the oscillator is 
measured to lie within a certain location in phase space, (a) 
LO = 41171.6 rad s~^. A pair of oscillation states emerges near 
ujhi- (b) uj — 41153.0 rad s"'^. As u decreeises, the two states 
move further apart in phase space, (c) uj = 41139.8 rad s~^. 
When u) < aj(,2, an addition attractor at the origin appears, 
(d) LO = 41124.8 rad s~^. With further decrease in oj, the 
occupation of the zero-amplitude state increases. 



values that differ by tt. When the oscillator is in the zero- 
amplitude state, there are large fluctuations of the phase 
as a function of time. The coexistence of three attractors 
in phase space is also illustrated in Figs. 13; and|31i for 
two other driving frequencies. 

We identify the residence time in each state before a 
transition occurs and plot the results as a histogram for 
one of the oscillation states (Fig. The exponen- 

tial dependence on the residence time indicates that the 
transitions are random and follow Poisson statistics as 
expected. From the exponential fit to the histograms, 
the transition rate out of each state is extracted. The 
transition rates out of the two oscillation states are mea- 
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FIG. 4: Histogram of the residence time in one of the oscil- 
lation states before switching occurs, at a driving frequency 
of 41163.0 rad s~^. The solid line is an exponential fit. (b) 
Logarithm of the transition rate as a function of inverse noise 
intensity. 



sured to be identical to within experimental uncertainty 
at all noise intensities. Figure plots the logarithm 
of the transition rate as a function of the inverse noise 
intensity. The transition rate varies exponentially with 
inverse noise intensity, demonstrating that the switch- 
ing is activated in nature. Using Eq. |21 we obtain the 
corresponding activation barriers at a particular driving 
frequency from the slope in Fig. ^p. Transitions out of 
the zero-amplitude state are also found to be activated 
and follow Poisson statistics in a similar manner. 

The above procedure is repeated to determine the ac- 
tivation barriers at other driving frequencies. Figure 5a 
shows the driving frequency dependence of the activa- 
tion barriers Ri and R2 for switching out of the oscil- 
lation states and the zero-amplitude state respectively. 
Ri is calculated using the average transition rates for 
the two oscillation states. At the high frequency end 
of Fig. only the zero-amplitude state is stable. As 
the frequency is decreased, two stable oscillation states 
(separated by an unstable state) emerge at ujbi. With 
increasing frequency detuning Acji — cjbi — w, the pair of 
oscillation states move further apart in phase space (Figs. 
3a and 3b) and Ri increases. At lo}j2, the zero-amplitude 
state becomes stable. The appearance of the stable zero- 
amplitude state is accompanied by the creation of two 
unstable states separating it in phase space from the two 
stable oscillation states. R2 initially increases with fre- 
quency detuning IS.UJ2 = ^b2 — a; in a fashion similar to 
Ri. Close to Wh2, Ri exceeds R2 and the occupation of 
the oscillation states is higher than the zero-amplitude 
state (Fig. As the frequency decreases, R2 con- 

tinues to increase monotonically while Ri remains ap- 
proximately constant. As a result, Ri and R2 cross each 
other at 41140 rad s~^, beyond which the occupation 
of the zero-amplitude state becomes higher than the os- 
cillation states. The dependence of the occupation on 
frequency detuning was also observed in parametrically 
driven atoms in magneto-optical traps 

In general, the activation barriers Ri and R2 depend 
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FIG. 5: Dependence of the activation barriers Ri (squares) 
and R2 (circles) on the frequency of parametric modulation, 
(b) Ri vs Aloi on a logarithmic scale. The lines are power law 
fits to different ranges of Acui . (c) R2 vs Au}2 on a logarithmic 
scale. 

on various parameters of the device, including the damp- 
ing constant, nonlinear coefficients, modulation ampli- 
tude and frequency. Nonetheless, at frequencies close 
to the bifurcation points, theoretical analysis indicates 
that the activation barriers exhibit universal scaling, with 
R1.2 = 1^^61,2 — For pitchfork bifurcations in a para- 
metric oscillator that involve merging of three states, ^ is 
predicted to be 2 |Q| . Figures 5b and 5c show the depen- 
dence of R1.2 on frequency detuning Awi_2 on logarith- 
mic scales. At small detuning, both activation barriers 
show power law dependence on detuning. The critical 
exponents are measured to be 2.0 ±0.1 and 2.00 ± 0.03 
for Ri and R2 respectively, consistent with theoretical 
predictions. This quadratic dependence of the activa- 
tion barrier on detuning near the bifurcation points is 
system-independent and is expected to occur in other 
parametrically-driven, nonequilibrium systems such as 
electrons in Penning traps [^f and atoms in magneto- 
optical traps 0. Away from the vicinity of the bi- 
furcation point, however, the variation of the activation 
barrier with frequency detuning is device-specific. Fig- 
ures 5b and 5c show crossovers from the quadratic depen- 
dence to different power law dependence with exponents 
1.43 ± 0.02 and 1.53 ± 0.02 for Ri and R2 respectively. 
These values are distinct from the exponents obtained in 
parametrically driven electrons in Penning traps be- 
cause the nonlinearity and damping are different for the 
two systems. 

Recent theoretical predictions indicate that the sym- 
metry in the occupation of the two oscillation states in a 
parametrically driven oscillator will be lifted when an ad- 
ditional small drive close to frequency aj/2 is applied [l^ . 



A number of phenomena, including strong dependence of 
the state populations on the amplitude of the small drive 
and fluctuation-enhanced frequency mixing, are expected 
to occur. Further experiments are warranted to test such 
predictions and reveal other fluctuation phenomena in 
parametrically driven oscillators. 

Parametric pumping is widely used to improve the 
sensitivity of micromechanical detectors by mechanically 
amplifying a signal ,18i] or by reducing the resonance 
linewidth in viscous environments |l9l |. The sharp jump 
in oscillation amplitude at the bifurcation points is uti- 
lized for accurate determination of device parameters 
the dissipation increases in a resonantly-driven 
Duffing oscillator, the natural resonance linewidth be- 
comes very broad and the hysteresis region shrinks for 
comparable oscillation amplitude. Parametrically driven 
oscillators, on the other hand, maintain the sharp jump in 
response at the subcritical bifurcation point even for large 
damping, ft is not necessary to increase the oscillation 
amplitude provided that the stronger parametric pump- 
ing compensates the energy loss due to damping 20l|. 
Therefore, parametric oscillators are particular useful for 
sensing in liquid or gaseous environments. Apart from 
the relevance to other parametrically driven nonequilib- 
rium systems 0) IE j the comprehensive study of the 
dependence of the transition rate on frequency reported 
here may prove useful in sensing applications with para- 
metrically driven micromechanical devices. 

We thank M. f. Dykman and D. Ryvkine for useful 
discussions. 
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